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Abstract 

We study the asymptotic Hmit of solutions to the barotropic Navier-Stokes sys- 
tem, when the Mach number is proportional to a small parameter e — )• and the 
fluid is confined to an exterior spatial domain il^ that may vary with e. As e — )• 0, 
it is shown that the fluid density becomes constant while the velocity converges to 
a solenoidal vector field satisfying the incompressible Navier-Stokes equations on a 
hmit domain. The velocities approach the limit strongly (a. a.) on any compact set, 
uniformly with respect to a certain class of domains. The proof is based on spec- 
tral analysis of the associated wave propagator (Neumann Laplacian) governing the 
motion of acoustic waves. 
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1 Introduction 

There is a vast number of mathematical models of incompressible fluids that can be iden- 
tified as a singular limit of more complex systems describing the motion of compressible 
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and/or heat conducting fluids, see the review papers by Gallagher [TS], Masmoudi 
[25], or the monograph [15] and the references therein. In many cases the resulting sys- 
tem is driven by an exterior force related to gravitation of rigid objects outside the fluid 
domain whereas, at the same time, the fluid occupies the whole physical space R^, see the 
study of the Oberbeck-Boussinesq approximation by Brandolese and Schonbek [3]. Such 
a situation, if physically relevant, can be viewed as a singular limit, where the primitive 
system is posed on a family of exterior domains fi^, with R^\Qi; being the rigid body(ies) 
acting on the fluid by their gravitation, and — )■ i?^ as e — )■ in a certain sense. Our goal 
in the present study is to develop a method for studying the incompressible limits, with 
the Mach number Ma = e — )■ 0, on a family of exterior domains fl^ C R'^ varying with 
e > 0. In particular, we identify a class of domains giving rise to uniform convergence 
of the fluid velocities, independent of the speciflc shape of their boundaries. To this end, 
we adapt the technique introduced in [H], based on spectral theory for the correspond- 
ing acoustic wave propagator - the Neumann Laplacian on Q^. The adaptation leans on 
delicate estimates of the associated Hemholtz projections based on the results by Farwig, 
Kozono, and Sohr in [12] and |13j . 

For the sake of simplicity, we focus only on the mechanical aspects of the fluid motion, 
ignoring completely the effect of temperature changes. Accordingly, we consider the 
compressible Navier-Stokes system in Eulerian reference coordinates: 

+ div^.(^u) = 0, (1.1) 
dtigu) + div^(^u ® u) + \v^p{g) = div^S(V^u), (1.2) 

S(V^u) = fi (v^u + V^u - ^div^ul^ , > 0, (1.3) 

where the unknown functions are the mass density g = g{t, x) and the vector velocity 
u = u(t, x), with t e (0,T), and x G fie- The symbol p = p{g) denotes the pressure - a 
given function of the density - whereas S stands for the viscous stress tensor. The small 
parameter e represents the Mach number tending to zero in the asymptotic limit. 

The system (11. 11 - 11. 3j) is supplemented by the (acoustically hard) complete slip bound- 
ary conditions 

u • nlen^ = 0, [S ■ n] x n\en^ = 0, (1.4) 

where the n denotes the (outer) normal vector to dfl;;- Moreover, as the fluid occupies an 
exterior domain, the behavior of g, u at inflnity must be specifled: 

g ^ g > 0, u — )■ as |a;| — )■ oo, (1.5) 
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where ^ is a constant. 

In the singular hmit e — 0, the fluid is driven to incompressihility as the speed of 
sound tends to infinity. The incompressible limit for system (11.11 - II. 3p as well as related 
problems have been studied by many authors, see Alazard Klainerman and Majda 
[23], Lions and Masmoudi |30], or Schochet |10], to name only a few. Similarly to Lions 
and Masmoudi [SU], [SI], our approach is based on the concept of weak solutions to the 
Navier-Stokes system, where the presence of viscosity plays a crucial role. As is well known 
(cf. [30]), convergence of solutions of the system (11.11 - [L3l) to the incompressible limit 
may be disturbed by the presence of rapidly oscillating acoustic waves, here represented 
by the gradient part of the velocity field. Since the physical domains Vt^ are unbounded, 
however, we expect that acoustic waves disperse leaving very fast any bounded part of the 
physical space as was observed by Alazard [T], Bresch and Metivier [3], Isozaki [2T], and 
[T7] . The main novelty of this paper is the fact that the physical domains f2e are allowed 
to change their shape together with the Mach number. In particular, we show that the 
rate of convergence is uniform within a certain class of domains specified in Section 12.21 
below. An interesting aspect of the problem is the boundary behavior of the limit velocity 
field U. As shown in P], the slip boundary conditions (ll.4p may give rise to the more 
standard no-slip condition 

U|of, = 0, (1.6) 

or to the kind of friction-driven boundary conditions identified in [7]. 

In order to see the principal difficulties involved, we rewrite the Navier-Stokes system 
in the form of Lighthill's acoustic analogy [2Z], |28] : 

edtr + div^V = in (0, T) x n,, (1.7) 

edtV + p{g)V^r = ediv^L in (0, T) x Q„ (1.8) 
supplemented with the boundary condition 

V-n|an, = 0, (1.9) 

where _ 

r = ^^, V = gu, (1.10) 

and L is the so-called Lighthill's tensor, 

L = §-^u®u - ^{p{g) - p'(:g){g -g)- pig))i. (1.11) 

Applying, formally, the Helmholtz projection to (11. Sp we obtain the wave (acoustic) equa- 
tion 

edtr + A^ = 0, (1.12) 
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edt<^ + p'{Q)r = eA-]vdiv^div^L, (1.13) 

V,$-n|an, =0, (1.14) 

for the acoustic potential $ = A^j^rdiv^V, where the symbol A^^n denotes the standard 
Laplace operator supplemented with homogeneous Neumann boundary conditions on dQ^- 
Accordingly, our main task consists in: 

• estimating the forcing term A~jydiVa;diVa;L as well as the initial data in terms of a 
suitable power of the operator A^^n, cf. |T7]; 

• evaluating local rate of decay of solutions to the wave equation f ll.l2p . fll.lSp . 

In general, bounds of Ajj^rdiVj-div^jL in terms of Ag tv depend on the shape of dfl^, where 
the latter must be at least of class C^'^ to recover the standard 1^^'''— theory, not available 
on less smooth, say, Lipschitz domains, see Grisvard [I9l . On the other hand, our method 
is applicable to families {fl^}e>o, whose smoothness parameters blow up for e — )■ 0. In 
particular, they may approach a less regular domain in the asymptotic limit and/or their 
boundaries may oscillate similarly to [S], see Section [2.21 As a result, the forcing term 
involving Lighthill's tensor become unbounded for e — 0, and this defect must be com- 
pensated by uniform dispersive estimates of order ^/e which we achieve in the spirit of a 
result due to Kato [23] . 

It is interesting to note that similar results on bounded domains, supplemented with 
the no-slip boundary condition (11.61) . where convergence of the velocities is enforced by 
a viscous boundary layer (see Desjardins et al. [Hj) seem very sensitive and much less 
stable with respect to domain perturbations, in particular they completely fail on balls. 

The organization of the paper is as follows. In Section [2|, we recall some known 
facts concerning the compressible Navier-Stokes system (ll.ll - [r!3|) . including the available 
existence theory, introduce the principal hypotheses concerning the admissible class of 
domains, and state our main result. Section [3] reviews the standard uniform bounds on 
solutions to (11.11 - II. 3p in (0, T) x fl^ independent of the singular parameter e — )■ 0. In 
Section HI we introduce the acoustic equation and identify the terms in Lighthill's ten- 
sor. In particular, we deduce estimates on Lighthill's tensor in terms of e, based on the 
ly^'P— theory for the Neumann Laplacian Ag ^v defined on fi^. Section [5] deals with the 
spectral theory of the operator — A^ tv in L'^i^e)- We introduce the associated spectral 
measures and show local decay of acoustic waves of rate ^/e. Finally, the incompress- 
ible limit is performed in Section [61 In addition, we shortly discuss the limit boundary 
conditions in the spirit of [7]. 
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2 Preliminaries and main result 

Throughout the paper, the pressure p is a continuously differentiable function of the 
density such that 

p e C[0, oo) n C^(0, oo), p{g) > for all ^ > 0, Jim = > 0, (2.1) 

for a certain 7 > 3/2. 

Multiplying the momentum equation (11. 2p by u and integrating by parts leads to the 
energy inequality: 

/ E,i0,u)iT,-) dx+ [ §(V,u) : V,uda;<0, TG (0,T), (2.2) 
at J^e Jn^ 



where we have set 
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EeiQ, u) = -g\u\' + - {P{q) - P\q){q -q)- Pm , 
2 s'^ 

P(g) ^ g r dz. (2.3) 

As P"{q) = p'{q)/q > 0, the function P is strictly convex and P{g) — P'{'q){q — '^) — 
P i'g) ^ c{g — 'g)'^ provided g ^'g. Consequently, with initial data of the form 

^'(0, •) = go,s = Q + £^^o,e, u(0, ■) = Uo,e, (2.4) 

with 

lko,e||L2{Q^) + lko,e||L°°(ae) + II Uo,e || L2(n^.R3) < C, (2.5) 

we get the total (mechanical) energy associated to the initial data 

Ee{go,s, Uo.e) = Q^'0,e|uo,£p + ^ {PiQo,e) " ^'(^)(^'0,£ ^ Q) - ^(^))) dx, 

bounded uniformly for e — )■ 0. 
2.1 Weak solutions 

We say that a pair of functions g, u represents a weak solution to the Navier-Stokes system 
(O-O]), with the boundary conditions (^^, (O]), and the initial data (|231) if: 
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g>0,ig-g)e {L'^ + L"<){Vl,) , u e L'^{f),T]W^^'^{n^]R^)) such that u ■ n|af^^ = 
the equation of continuity (11.11) is satisfied in the sense of renormahzed solutions: 

/7 



(2.6) 



dxdt = - {go^e + b{0o,e))^iO,-) dx, 

for any test function ^ G C^{[0,T) x H^), and any b E C~[0, 00), 6' G C^[0, 00) 
the momentum equation (II. 2p holds in the sense of the integral identity 



T r r I 

gu ■ dt(f + gu ig) u : V^ip + —p{g)diw^ip 



Jn^ 

T 



dx dt 



(2.7) 



= / / S(V^u) : V^ip dx dt - / ^io,£Uo,£ ■ V5(0, ■) dx, 
for any test function (f G C^{[0,T) x Q^; R^), V x'^ ■ ^\d^e = 0? 



the energy inequality 
/o 



- r E,{g,u)dt^ dt + r ^ I §(V^.u) : V^u dx dt 
Jo Jo Jn^ 
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< 



^e(^'0,e,Uo,e)V^(0) dx. 



holds for any G C~[0, T), t/- > 0. 
2.2 Admissible domains 

Motivated by |8] , we introduce a class of admissible domains that allows for "oscillating" 
boundaries. As first observed by Casado-Diaz, Fernandez-Cara, and Simon [9], such 
a family of domains may give rise to the no-slip boundary condition (II. 6p for the limit 
velocity field. A general description of all possible limit boundary conditions was obtained 
in 

Introducing a cone 



C(x, uj , 5, ^) = {y E \ Q < \y — x\ < 5, {y — x) ■ ^ > cos((X')|y — x|}. 
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with vertex at x, aperture 2a; < vr, height S, and orientation given by a unit vector ^, we 
say that satisfies the uniform 5— cone condition if for any xq G Sfi^, there exists a unit 
vector ^ such that 

C{x, u, S, ^xo) C ^£ whenever x G fie, \x — Xo\ < 5, 

see Henrot and Pierre [201 Definition 2.4.1]. 

Assume we are given a family of domains {Qs}e>o complying with the following hy- 
potheses: 

• Qe G is an exterior domain with C^— boundary for each fixed e > 0; 

• there is a d > such that 

R^\n, C Bd = {x e \ \x\ < d} for all e > 0; 

• satisfy the uniform 5— cone condition with 6 > (and u) independent of e; 

• for each xq G dQe, there are two (open) balls -Br[a^i] = {x;\x — Xi\ < r} C fi^, 
Br [xe] G R^ \fle of radius r > CbS^ such that 



Br[x.-]nBr[Xe]=Xo, (2.9) 

with Cfe > 0, /3 > independent of e. 
The above hypotheses give rise to the following properties enjoyed by the family 

{^e}e>0: 

• Uniform extension property (see Jones f2^). 
There exists an extension operator E^, 

E, : W''P{Q,) -> W''P{R^), E,[v]\n. = II^.H Iki..(ij3) < c\\v\\w^.,^n^), (2.10) 
where the constant c is independent of e — 0. 

• Uniform Korn's inequality (see 16, Proposition 4-^])- 

Let V G W^''^{Sl, n B] R^), and M C n 5 such that \M\ > m > 0, where 5 is a 
bounded ball. Then 

l|v||^i,2(n,ni?;i?3) < c(m) ( V^v + V*.v - ^div^vl + j |v|2 dx ) , 

(2.11) 

with c(m) independent of e — )■ 0. 
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• Compactness (see Henrot and Pierre 12^ Theorem 2.4-10]). 

There exists an exterior domain fi, satisfying the uniform 5— cone condition, and a 
suitable subsequence of e' s (not relabeled) such that 

1^^ + ^ as £ ^ 0. (2.12) 

For each xq G dVt, there is x^^o ^ d^e such that xq, in particular, 

R^\VLcBd. (2.13) 

For any compact K (zVL., there exists s{K) such that 

K ^VLe for all e < e{K). (2.14) 

Property (12.121) is important when studying stability of the spectral properties of the 
Neumann Laplacian A^jv with respect to e, see Arrieta and Krejciffk p]. Note that the 
limit domain need not be of class but merely Lipschitz, see Henrot and Pierre [20t 
Theorem 2.4.7]. 

2.3 Main result 

Before stating our main result, we introduce the limit problem - the incompressible Navier- 
Stokes system - satisfied by the limit velocity field U: 

div,.U = 0, (2.15) 
^((9tU + div^(U®U)) + V^n = /iAU, (2.16) 
in (0,T) X with the condition at infinity 

|U| ^ as |x| cx), (2.17) 

and the initial condition 

U(0,-) = Uo. (2.18) 

In the weak formulation, the decay condition f l2.17p is replaced by a single stipula- 
tion U G L^(0, T; W^^'^(fi, i?^)), the incompressibility constraint ( I2.15P is satisfied a. a. in 
(0,T) X f2, while the momentum equation fl2.16p . together with fl2.18p . are replaced by a 
family of integral identities; 

r I {-q\J ■dtV + -Q'\J ®\J : V^<^) dx dt (2.19) 
Jo Jn 
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= / / Va;U : VxV dx dt - / ^Uo ■ V5(0, ■) dx, 
JO Jn Jn 

for any test function (p G C^([0,T) x fi; i?^) satisfying diVx^ = 0. 

Note that we have dehberately omitted to specify any boundary conditions on dQ. 
In the low Mach number hmit, we can easily show that the impermeability condition 

■ n|af7^ = gives rise to the same property U • n\sn = for the limit velocity field 
whereas the boundary behavior of the tangential component of U may be quite complex 
depending sensitively on the asymptotic shape of the boundaries dQs, cf. |2]. Sufficient 
conditions for U to satisfy the no-slip condition 



U\an = (2.20) 



will be discussed in Section 16.21 
Our main result reads as follows. 



Theorem 2.1 Suppose that a family of domains {^^e}e>o belongs to the class specified in 
Section IKR with 

o<P<l, 

where is the radius of the halls in Ii2.9\) . Let {f)e,U£}e>o be a family weak solutions to 
the compressible Navier-Stokes system U.l\ - [773]) in (0,T) x Q^, supplemented with the 
initial conditions \2.Ji\), where 



ro,e = ^ -> ^0 weakly in L'^{R^), \\ro,e\\L^{R3) < c, > 0, 

Uo,e — Uo weakly in L^{B?] R^), 

and the pressure satisfies h2.1\) with 7 > 3/2. 
Then 

ess sup II •)- ^11 (L2+i5)(n^) ^ as £ /or 1 < g < min{7, 2}, (2.21) 

iG(o,r) 

||u£||L2(o,r;vyi'2(Q^;R3)) < c, (2.22) 
and, at least for a suitable subsequence, 

Ue -> U zn L2((0, T) X K] R^) for any compact K C fi^, (2.23) 

where 

U G L2(0, T; W^'\^- R^)), U ■ n\dn = 0, 
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is a weak solution of the incompressible Navier-Stokes system l{2.15\ - 2.11} in (0,T) x Vt, 
emanating from the initial data 

Uo = H[uo], 



where VL is the limit domain identified through l{2.12\ - \2.14\ ) and H denotes the standard 
Helmholtz projection in Q. 

The rest of the paper is devoted to the proof of Theorem I2.1[ Note that existence of 
solution to the compressible Navier-Stokes system, in the framework of weak solutions, was 
first established in the seminal work of Lions [29] for 7 > 9/5, and the result then extended 
in [16] to the "technically" optimal range 7 > 3/2. A detailed discussion of various choices 
of boundary conditions, including the case of exterior domains, may be also found in the 
monograph by Novotny and Straskraba [3^. The problem of limit boundary conditions, 
here discussed in Section [HI was studied by Casado-Diaz, Fernandez-Cara, and Simon 
[S] in the periodic setting, and later extended in [7], [S]. Finally, we remark there is an 
analogue of Theorem 12 .ll for a family of bounded domains obtained by completely different 
methods (see [6]), based on analysis of a viscous boundary layer similar to Desjardins et 

al. [n]. 



3 Uniform bounds 

All uniform bounds presented below may be viewed as a direct consequence of the energy 
inequality fl2.8p . To begin, similarly to [IS], we introduce the essential and residual part 
of a function as 

[^e]ess X(^?e)^e) [^e]res h /^ess, 

where 

X e Cr^(0, 00), 0<x<l, x = linan open neighborhood of 'g. 
3.1 Energy bounds 

Since the initial data satisfy (12 ■4p . (12.51) . the initial energy i?£(f)o,£, uo^^) in (12. 8p remains 
bounded uniformly for e — )■ 0, where we have used hypothesis (12.10 . Consequently, we 
deduce the following list of estimates: 

ess sup / f?e|u£p dx < c, (3.1) 

ie(0,T) 
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ess sup 



and 



/ 


'Qe - Q' 




e 



dx < C, ess sup / [QeVres ^-"^ — ^"^^i / "'-res dx < E^C, 

(3.2) 



T 

JD,e 



dx dt < c, 



(3.3) 



where the constants are independent of e — )■ 0. 
Moreover, relation (13. 2p immediately yields 



ess sup 

ie(o,T) 





'Qe - Q' 






e 


res 



2-q 



< ce 1 for any 1 < q < min{7, 2}, 



(3.4) 



which, together with (13. 2p gives rise to (I2.2ip . 

Finally, since the family of domains {^£}e>o admits the uniform Korn's inequality 
(12. lip , we can combine (13. 2p . (13. 3p to conclude that 



|Ue||L2(0,T;iyi>2(Oe;i?3)) < C, 



(3.5) 



uniformly for e — )■ 0. 



3.2 Convergence 

As the family {Qs}e>o possesses the uniform extension property (I2.10p . we may assume 
that 

Ue^V weakly in ^^(0, T; W^^^R^; R^)). (3.6) 

Moreover, by virtue of the uniform bounds (13. 2p . (13. 5p . we can perform the limit in 
(12. 6 p to obtain 

div^U = a.a. in (0,T) x fi, 

and, similarly, one can pass to the limit in the weak formulation of momentum equation 
(12.70 to deduce that 

/ / (■^U ■ dfcp + gu^u : Vx^) dx dt 
Jo Jn 

= /i / / Vx'U : Vx^f dx dt - / ^Uo ■ (^{0, ■) dx, 
Jo Jn Jq 

for any test function ip G C^([0,T) x Q; R^) satisfying divxf = 0, where the symbol 
® u denotes a weak limit of {fPeU^ (g) U£}£>o- Accordingly, in order to finish the proof 
of Theorem 12.11 we have to show 



gu0 u = gU U, 
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or, equivalently, the strong convergence of the velocities claimed in fl2.23p . This will be 
our goal in the remaining part of the paper. We remark that, by virtue of (12.121 - I2.14p . 
it is easy to check that the the limit velocity field satisfies the impermeability condition 

U • n\on = 0, 

in a weak sense. 

4 The acoustic equation 



In this section, we introduce a weak formulation of the acoustic equation (I1.12I - [LT^ and 
discuss its basic properties. 

4.1 Weak formulation 

With 

re = , = QeUe, 

the Navier-Stokes system (12. 6p . (12. 7p can be written in the form: 

r [ (eredtif + Ve ■ V^v?) dx dt = - [ eroM^^ ■) dx, (4.1) 
for any ip G C~([0,T) x U), 

[ [ (eVe ■ dt^p + p'{g)rsdiv^ip) dx dt = - [ e^io.eUo,^ ■ (^{0, ■) dx (4.2) 

L In (^^*^^^"^) ■ '^^'^ ~ ^^^^^ ® ■ '^^^) dt 

~ lo In ^ (^(p(^'e) - -q)- P{p))div^(p^ dx dt, 

for any cp G ([0, T) x H,; R^),^- n\an. = 0. 

Now, thanks to the slip boundary condition (II. 4p . observe that A j]^ where A^^n 
is the Neumann Laplacian in , is an admissible test function in (14. 2 p . Consequently, we 
obtain 

r I (e^edtV - p'i-gyev) dxdt = - [ £f?o,eUo,. ■ V.A;]^[(^(0, ■)] dx, (4.3) 

12 



- r f (£§(V.u,) : VlA-'M - ^QeUe ® u, : V^A-jvb]) dx dt 

for all (p G C^([0,T) x where $e is the acoustic potential, meaning, 

V, = H,[V,]+V,.<I>„ 

where denotes the standard Helmholtz projection in fi^. Note that for ip G C^([0, T) x 
He) and e > fixed, the test function Va;A~jv[(y9] is continuously differentiable in il^ and 
belongs to the space C~([0,T); W^'^in^; R^))- 

The Helmholtz projection H£[v] of a function v G L^^Qe', R^) is defined as 

H,[v] = V- V.'l', 

where $ G ^^'^(ne) is the unique solution of the problem 

/ V:,$ ■ V^ip dx= [ V ■ V^ip dx for all ip G 

in other words, at least formally, 

A$ = div^v in G fi^, Va;$ ■ n\Qn, = v ■ n\o^^, |$| for |x| oo. 

Here, the symbol D^'^{Qe) denotes the homogeneous Sobolev space - a completion of 
C^{fls) with respect to the norm || Va:$||LP(Cj;i?3)- We have Sobolev's inequality 

IML^n,) < c(p)||V^.$|Up(n,;ij3), q = j—^ for any 1 < p < 3, (4.4) 

for any $ G where, since {^^£}e>o admits the uniform extension property fl2.10p . 

the constant c(p) is independent of e. 
Finally, the equation (14. ip reads 

r f (eredtip + V^.$, ■ V^y?) dx dt = - [ ero^e^{0, ■) da; (4.5) 
Jo Jn^ ^ ' J^e 

for any ip G C^([0,T) x fi). Equations (14. 3p . (14. 5 p represent a weak formulation of the 
acoustic equation (I1.12p . (I1.13p . with the Neumann boundary condition (I1.14p implicitly 
included through the class of test functions in (14. 5p . 
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4.2 Uniform bounds, part I 

It follows from the uniform bounds established in f l3.ip . (I3.2p . (13 ■4p that 
satisfies 

ess sup ||[r4,J|^2(f^ ) < c, (4.6) 

and 

ess sup ||K]i.csllL9m ) - ^ " ^ 1 < 5' < min{7,2}. (4.7) 

te(o,T) ^ 

Similarly, 

where, in accordance with (13. ip . (13. 2p . 

ess sup ||[Ve]css|U2(Q^.R3) = ess sup \\[^/g^]ess^/0^Ue\\L2{n,■,m) < c, (4.8) 
tG{0,T) i6{0,T) 

and 

ess sup ||[Ve]res|U9{C,;R3) = ess SUp \\[y/g^]rcs\^U^\\ Li{Q,;R3) < Ce'^^'^ , q = —■ (4.9) 

ie(o,T) te(o,T) 7 + i 

It remains to find suitable bounds on the forcing terms in acoustic equation (14.31 - 14. Sp . 
To this end, we need the elliptic estimates for the Neumann Laplacian a? discussed in 
the next section. 

4.3 Elliptic estimates and Helmholtz decomposition in 

In order to control the forcing terms as well as the initial data in the acoustic equation, 
we need bounds on V^v in terms of A^^Tvf'f]- As the curvature, represented by the radius 
of the balls in (12. 9p . is not uniformly bounded, the ly^'''— elliptic bounds may "blow-up" 
for e 0. 

4.3.1 iy2,p_bo^jjjjg 

In order to obtain W^^'^— bounds, we consider the rescaled family of domains 

a = ^^e, (4.10) 

14 



where the exponent /3 > is the same as in fl2.9p . Accordingly, the rescaled domains 
are of uniform C^-class, in particular, the standard elliptic theory yields 

ll^'^llLP(n,;i?3x3) < c{p) {\\A.,v\\^^^-^^ + for any 1< p < oo, (4.11) 

for any v G C^{fls) satisfying V^f ■ n|gj^ = 0. It is important to notice that, by virtue of 
the hypotheses introduced in Section I2.2[ the constant c{p) depends only on the rescaled 
radius of the balls appearing in fl2.9p . 

Consequently, returning to the original domains fig we may infer that 

||V^t;||LP(Q^;R3x3) < c{p) (^\\A^v\\LP{n,) + ^||^^||LP(f7,)) for 1 < p < oo, (4.12) 

for any v G C^^fls) such that VxV ■ nlan^ = 0, with c{p) independent of £ — 0. 

4.3.2 Helmholtz decomposition 

Consider the family of rescaled domains Qg introduced in fl4.10p . with the associated 
Helmholtz projections H^. By virtue of the result by Farwig, Kozono, and Sohr [12j, we 
have 

\\^eH\{LpnL^m,R^) ^ <P)M{LvnL2m,R^) any 2 < p < oo, (4.13) 

where, similarly to the preceding part, the constant c{p) depends only on Cb- Going back 
to the original domain we therefore obtain 

||He[v]||(LpnL2)(f7„ii3) < £'^^^"1^^ c{p)\\\\\i^LpnL^){n„R3) for any 2 < p < oo, (4.14) 

uniformly for e ^ 0. 

Similarly, by means of a duality argument, 

\\iie[^]\\iLP+L^)in„R3) < e'^^"^'^^ c{p)\\\\\(^LP+L^){n,,R^) for any 1 < p < 2. (4.15) 

The estimates 04.140 . (14.150 arise also in problems of homogenization and need not be 
optimal, cf. Masmoudi [33], [34j. On the other hand, as the limit domain Q may be only 
Lipschitz, it is not surprising that the L^— bounds in (14.140 . (I4.15P blow up for e — )• 0. 

4.4 Uniform bounds, part II 

With the bounds established in the previous section at hand, we are able to control the 
forcing terms in the acoustic equation (14. 3p . To begin, relation (14.120 implies that 

/ §(V,u,) : V^.A-jvM dx 
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therefore, by means of the Riesz representation theorem, 

r [ S(V.u,) : V^.A,"jvM dxdt = ^ r [ (fIv + F,2(-A,,A.)-iM) dx dt, (4.16) 
where 

l|i^£ilL2((o,T)xf7,) < c, 2 = 1,2, uniformly for e ^ 0. (4.17) 
Similarly, we can write 

g^u^ (g) Ug : VIA^)^[lp] dx 

= / [^leJessUe O : V^A;;]v[v9] dx + / [^/g^]res^/Q^U, (g) : V^,A;;jvb] dx, 

where, by virtue of fl4.12p . 

/ [&]cssU£ ® : v^A;;jvM dx 

< ||[&]essU£||L2(Q^^^3)||Ue||i6(n^./j3)||V^.A;^]v[V2]||L3{Q,;_R3x3) 
< c||[^£]essUe||L2(Q^,fl3)||Ue||L6(o,;i?3) (^||</?||L3(f^^) + || (- A^,^)"^ || i3(n,)^ , 

and, by interpolation, the uniform extension property, and Sobolev's inequality, 

||¥'IU3(f^^) < ci i\\ip\\L^n,) + WvWh^n, 



< C2 + ||V^(/?|U2(f^^)j = C2 + ||(-A,,Ar)l/2[^]||^^^^^ 

while, by the same token, 

\\{-Ae,NVM\\L^in.) < ci (||(-A.,A.)-Mv^]||L^(n.) + II ("A.,^)"^/' IU^(n. 
Consequently, there exist functions Fl, i = 3, . . . , 6, 

||^el|L2((o,r)xn,) < c, i = 3, . . . , 6 uniformly for e ^ 0, (4.18) 

such that ^ 

/ / [&]essu, ® u, : V^A;^[yp] dx dt (4.19) 
JO Jn^ 
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Furthermore, in accordance with (14. 9p . 

[&]rcsU£ ® : VlA-]^[(p] dx < £^/^c||u||i6(n,;fl3)||V^A-jv[y5]||L'-(f7,;i?3x3), 

where, by virtue of f l4.12p . 



with 



-'+1 ' 1+1=1. 



,-1 



27 6 r 

In what follows, we suppose 7 < 2, r > 3 as, otherwise, the estimates are the same as 
in fl4.19p . Thus, applying once more fl4.12p . we obtain 



and, similarly, 



r2/3 



-A 



e,N) 



Consequently, we get the same result as f l4.19p provided max{-, 2} — 2/^, in particular if 
(3 < 1/4. 

Thus, we may infer that 



/ / QsUs ® Ue : V^A^ jv[v9] dx dt 
Jo J^le 



(4.20) 



1 

e^^ Jo Jn, 



+FH-A,.Nm+F'(-A 



dx dt, 



where satisfy (I4.18p . Note that the same symbol F^ may stand for different functions 
than above. 
Finally, 
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< 



Li(Q.) 



where, in accordance with the uniform extension property, 

||v'IU°°(n,) < c(\\'Vx^\\L(^{n,,R^) + ||</'||l6(q,)) < c (\\'Vlip\\L'2{n,,R3x-i) + W^x^Wi^in^-Ra) 
and, by virtue of fl4.12p . 



Seeing that, by virtue of the uniform bounds estabhshed <\3.2h . (13.40 . 



ess sup 

te(o,T) 



we conclude that 



1 

7^ 



where 



^£ IU2((o,r)xo,) < c, i = 8, 9, 10 uniformly for e -> 



(4.21) 



dx dt. 



(4.22) 



4.5 The acoustic equation revisited 



Using relations (14.81) . (14. 9p . together with (I4.15p . we can write the acoustic potential $£ 
in the form 



where 



ess sup ||$]||Dl.2(n,;/?3) < c, 

te(o,T) 

ess sup ||$^||2,2(f7^.^3) — )■ as £ — )■ 0. 
te(o,r) 



(4.23) 
(4.24) 
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In view of the uniform bounds obtained in the previous section, the acoustic equation 
( I4.3p . (14. 5p can be written in the concise form 

r [ (er.dtif + V^$, ■ V^<^) dxdt = - [ eroM^^ O (4.25) 

for any G C~([0,T) x H), and 

f [ (e^A^ - P'{0)re^) dxdt = - [ e(^oMO, dx (4.26) 

Jo Jn, ^ 

+Gti-A,,r,M + Gli-A,,Nr'M)dx dt, 
for all If G C^{[0,T) x U), V^if ■ n\dQ^ = 0, where 

ML^m + ||(-A,,^)-^/'[$o,e]||L^(c.) < c, (4.27) 
\\Gl\\mn^)<c, ^ = 1,...,5, (4.28) 

uniformly for e — )■ 0. 

5 Spectral analysis of Neumann Laplacian on varying 
domains 

As observed, the Neumann Laplacian jv plays a crucial role in the analysis of acoustic 
waves. We recall that — A^^jv niay be viewed as a non-negative self-adjoint operator on 
the space L'^^fl^), with 

P(-A,,jv) = {we W'^^Qe) I ^xw ■ V,(^ dx = / dx 
I Jn^ Jn. 

for all if G C^lQe) and a certain g G L^(fie)| , —A^^nw = g. 
Since the boundaries dfl^ are regular, the standard elliptic theory yields 

I?(-A,,^) = {we W'^^in,) I V.w ■ n\an. = o} . 

We denote by {Vs,x}\>o the associated family of spectral projections. The following 
analysis is a slight modification of [HI Section 2] , similar problems were studied by Ranch 
and Taylor 
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5.1 Spectral measures 

Our goal is to obtain dispersive estimates, and, in particular, local decay of acoustic waves, 
with a rate independent of the scaling parameter e. To this end, we introduce the spectral 
measure fie,<fi associated to a function ip G L'^{^e) through Stone's formula (see Reed and 
Simon |39i Theorem VII. 13]) 

IJ,e,^ia,b) (5.1) 
= lim lim / /( — 3— ^ | \(p];ip) dA, 

where the symbol (■; ■)^^ denotes the standard (complex) scalar product in L^(f2e). 

Now, the crucial observation is that we may perform the limit t] — )• 0+ in (15. ip as soon 
as V9 G C^(f2e) since the operators — satisfy the limiting absorption principle, see 
Leis Vamberg HH Chapter VIII]: 

Limiting absorption principle (LAP) 

The operators (1 + |xp)~*/^ o (— A^^jv — A ± ir])^^ o (1 + |xp)~'^/^ are bounded on L^(fie) 
for any s > 1 uniformly for A belonging to compact subsets of (0, oo) and rj > 0. 

Consequently, the spectral measure /i^^^, ip G C^{Qe), is absolutely continuous with 
respect to the Lebesgue measure on (0, oo) and Stone's formula (15. ip reduces to 



/^e,^(a, b) = {{w,^^ - w^^^^ ; <p)^^ dA, < a < 6, (5.2) 

where wf^ are solutions of the Neumann problem 

AwJ^ + Xwf:^ = (fin n^, V^wf^ ■ n|an, = 0, (5.3) 
uniquely determined by Sommerfeld's radiation condition 

Jiin r (j!)r ± Ia/a) wf^ = 0, r = \x\, (5.4) 

see Vamberg |!41j Chapter VIII]. 

Our goal is a uniform bound on the norm of the functions wf^ independent of the 
scaling parameter e. To this end, we fix i? > 2d so that dfl^ C -Br/2 for all e. Following 
we recall the explicit formula for the exterior Dirichlet problem 

+ A<A = OmR'\ Br, vf^^laB^ = vf^^, 
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supplemented with the radiation condition fl5.4p . namely 



txi^) = E E <Yri(^^ <P) /m ' for all xeR'\ Bn, (5.5) 

l=Om=-l h\ '{±VXR) 



where (r, 6, (f)) are polar coordinates, Y/" are spherical harmonics of order /, /i^'' are 
spherical Bessel functions, and 



oo / 



vU^) = E E for |x| = R, 

1=0 m=-l 

see Chandler- Wilde and Monk [lU], Nedelec [35] . 
Assume that G in fl5.3p is such that 

supp[v?] C Br. 

We claim that the functions wf^ solving f l5.3p . f l5.4p admit a uniform bound 

\\^fJ\LHB,nnn,] < cyh^in,), (5-6) 

with c independent of e — )■ 0, provided A belongs to a compact subinterval of (0, oo). In 
order to see (15. 6p . we argue by contradiction assuming the existence of sequences 



ifs e C^{Br n a), \\fe\\LHBi,nn,) = 1, 

Ae ^ A G (0,CX)), 

such that the corresponding (unique) solutions wf)^^ of (15. Sp . (15. 4 p satisfy 

\\wf^xJ\L^iB,nnn,) ^ oo for £ -> 0. 

Setting 



£ A — 1 1 di 1 1 

\\Ws,xJ\L^{B3Rnn,) 



we check that is the unique solution of the problem 

^^Ja + Kvf^x = II ± 11^" in ^s, V^vf^ ■ n\9u^ = 0, 

\\w^^xj\L^iB3Rnn,) 

satisfying 

\imr(dr±iJK] vf. =0, 
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and 

W^txjLHBsnnn,) = (5.7) 

Since 

A<A. = -A.<A, in B,R\Bn, 
the standard elliptic estimates yield 

vfx^ — )• in C^{K) for any m > and any compact K C -Bsij \ -Br, 

in particular, 

"^^A^ -)■ f^l" in C"^(9i?2i?) for any m > 0, 
which, together with formula (15. 5 p yields 

v^;^^ — )■ in C"^{K) for any m > and any compact K C \ -Br, 

where satisfy 

Au^ + Af ^ = in i?^ \ Br, Jiim r (^^ ± iVX) = 0. (5.8) 



Finally, we claim that satisfies 



A 

Av^ + At;^ = in n Bsr, V,v^ ■ n|ao = 0. (5.9) 
In order to see (15.91) observe that 

/ f V.<A. • V.^ - A,<,^V^ + ^ ^'"^ ) dx = (5.10) 

for any ijj G C^(f2 fl B^r). Since f^;^^ are bounded by (15.71) . we conclude that 

Wvf^xJw^'HBsRnQ,) < c, (5.11) 

uniformly for e — )■ 0. 

As fie possess the uniform extension property we may assume that (15. lip holds in R^, 
and, by virtue of (I2.12p . we may pass to the limit in (I5.10p for each fixed to obtain the 
desired conclusion (15. 9p . Moreover, by the same token, relation (15. 7p implies that 

WvxWLHBsnnn) = 1- (5.12) 
However, relations (15. 8p . (15. 9p imply f ^ = in contrast with (I5.12p . Thus we have shown 

(I5SD. 
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Summing up the previous discussion we may infer that 

< {(w-, - <,) ; < c(a, 6, ^) for any < a < 6, G (a), (5.13) 



uniformly for e 0. Moreover, repeating the arguments of the proof of fl5.6p . we conclude 
that 

}m({w-, - ; ^)^^ = ({w^ - wt) ; for any A > 0, ^ G C^{n), (5.14) 
where 

Aw^ + Xw^ = (fi inQ, V^w^ ■ n\dn = 0, 
lim r (dr ± wf = 0, r = \x\. 

Seeing that 

||V^(y9||L2(Q^.^3) = ||y^-Ae,Ar[v9]||i2(Q^.R3), (5.15) 

and _ 

(G(-A,,^)[^];^)^^^= / G(A)d/x,,^ 



for any if G C^{Qs), we deduce from f l5.13p . f l5.14p that 

||G(-A,,;v)MlkM(n.)<c(G,¥.), (5.16) 

and 

lf,,G(-A,,^)[V9] ^ lnG(-A^)M in ^^(7?=^), (5.17) 
for any G G C^{0, oo), G C,°°(l^). 

5.2 Uniform decay for |x| — > oo 

Our goal in this section is to establish the following decay estimate: 
Lemma 5.1 For G G C'^{Q, oo), ip G G'^{VL), supp[v9] C Br, we have 

I \x\^'\G{J^K^)[^f dx < c{G,s,R)M\l,^^) for any s > 0, 

J\x\>R * ^ ' 

uniformly for e ^ 0. 
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Proof: Extending G as an even function on it! we have 

G{^-^,,nM = ^ /^<5(^) (exp(i^-A,,jvi) + exp(-i^-A,,jvi)) b] di, 



where G stands for the Fourier transform of G. 
Denoting 



\s,R 



/ kl^kP'' dx, 

J\x\>R 



we get 



where 



< \ ||(exp(i^-A,,^t) +exp(-i^-A,,^t 

(exp(iy^A^i) + exp(-i^-A£,jv^)) W\ 
- R) \x\^' I (exp(i^-A,,^t) + exp(-iy'-A,,Aft 



s,R 



dt, 



2 



sgn {\x\ 



On the other hand, the wave operator 



dx. 



(^exp(iY'-A£,7vt) +exp(-iY'-A<.,7v^)j = 2 cos(^-A£,7v^) 
admits a finite speed of propagation 1, specifically, 

supp [(cxp(i^-A£,Art) +exp(-i^-Ae,Art)) [9?]] C BRm\] 

whence 

/ sgn+(|a;| -i?)|a;|2^ |(exp(i^-A,,wt) +exp(-i^-A,,^t)) [if] 
< + f |(exp(i^-A,,^rt) + exp(-i^-A,,^t)) [ip] 



dx 



dx 



^ {\t\ + Rn<p\\l.^a.y 
As G e C^{0, 00), we have {\t\ + RyC e L^{R) for any s, which completes the proof. 

Q.E.D. 
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5.3 Functional calculus - decay estimates 

Our ultimate goal in this section is the following result. 

Lemma 5.2 We have 



-t 



dt<£c(^,G)||vI/||i, 



for any if G C^{n), ^ G L\n,), and any G G C^{0, oo). 

Proof: We adapt the arguments of pJ:J. By virtue spectral theorem (see Reed and 
Simon [39], Chapter VIII]) we have 

(|exp (^iy^-A,,^^) [^],G(-A,,^)¥P^^ 

= l^'exp (^iv^^) G'(A)^,(A) d/i,,^(A), 
where fie,ifi is the spectral measure associated to the function (p, and 

Following Last we deduce 



^exp (i^-A,,N- ) [^],G'(-A,,^)[^] 



dt 



oo poo 



< eT0F / / exp (- ^''"^^^ G'(x)G'(y) d/i,,^(a;) d/x,,^(i/); 



JO 



therefore, by Cauchy-Schwartz inequality. 



-t 



[exp [i^-A.^N-j [%G{-A,,nM^ ^ 



dt 



POO I POO 

<eTV^J^ \^,{x)nj^ exp 



X 



Furthermore, 



exp 



X 



(5.18) 
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, ^, , exp dfi^Av) 

^0Jen<\^-^\<e{n+l) \ 6^ 4 J 



< sup / ldfle,^{y) 

n>0 ■Jen<\^-^\<e(n+l) „_r. 



I en<\^-V^\<e{n+l) \ 4 / 

Since only the points x belonging to supp[G'] are relevant in f lS.lSp . the length of the 
intervals 

In{x) = {y G [0, oo) I en < 1^ - v^l < e{n + 1)} 

never exceeds e, 

\In{x)\ < CgS. 

Thus we conclude combining f lS.lSp . with the uniform bounds on the spectral measures 
(see (I52D) established in fl5J[3|) . 

Q.E.D. 



6 Dispersive estimates, local decay of acoustic waves 

Returning to the acoustic equation fl4.25p . f l4.26p we can use the dispersive estimate es- 
tablished in Lemma 15.21 to show that 

|t ^ ^ <l>,G(-A,,jv)M dxj ^ in L\Q,T) (6.1) 

for any fixed G G C^(0,oo), ip G C^iVt). Indeed, by means of Duhamel's formula, we 
have 




+^'^^\ jl (exp {i^-^e,N^-^^ + exp (^-i^-A,,^^)) [H,{s)] ds, 

with 

He = Gl + (-A,,^)-i/2[G?] + (-A,,^)V2[g3] + (-A,,^)[G^] + {-^,,^)-'[G% 

(see fl4.26p ). where we have assumed, for the sake of simplicity, that p'{^) = 1. Conse- 
quently, the desired conclusion (16. ip follows from Lemma [5.21 as soon as /3 < 1/4, see |T3] 
for details. 
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6.1 Compactness in time of the momenta 

In view of estimate fl3.5p . the desired strong convergence claimed in (I2.23P follows as soon 
as we show that 

|t ^ ^Ve(-,t) ■ w da;| ^ |t ^^^U(-,t) ■ w da;| in L'^{0,T) for w G C^{n;R^). 

(6.2) 

Indeed relation (16. 2 p with (13. 6p imply that 

/ / f?e|u£p dx dt — 7- ^ / / |Up dx dt for any compact K C il, 
Jo Jk Jo Jk 

yielding (Km . 

In order to show (16.20 . we use Helmholtz decomposition to obtain 

/ ■ w dx = / He[^£Ue] -wdx- / $ediVa;W da; 
Jn^ Jn^ Jn^ 

= / ^eUe ■ He[w] dx - $ediv^.w dx 

Jn^: J^e 

= / QeMe ■ H[w] dx + / ^^u^ ■ (He[w] - H[w]) dx - / $ediv^.w dx, 

jQe J^E J^e 

where, in accordance with (14.20 and the standard Aubin-Lions argument, 

|t ^ ^ ^^u^ ■ H[w] dxj ^ |t ^ U • w dxj in ^^(0, T). (6.3) 

Here, we have extended H[w] and H£[w] by zero outside Q and il^, respectively. 
Furthermore, we write 

/ f5,u,-(H,[w]-H[w]) dx 
Jn^ 

= I (& - ^)u. ■ (H,[w] - H[w]) dx + ^ / u, ■ (H,[w] - H[w]) dx. 

In view of estimates (14. 6p . (14. 7p . and boundedness of H^, H in (see (14.141) . (14.151) ). we 

get 

\t ^ {Qe - -Q)^e " (H^M - H[w]) dxj ^ in L2(0, T). (6.4) 

Moreover, as 

He[w] ^ H[w] weakly in L^iVt] R^), 
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we get, by virtue of f l2.22p . 

[t^ ■ (H,[w] - H[w]) dxj ^ in ^^(0, T). (6.5) 

Finally, 

/ $£div^wdx=/ $eG(-Ae,^)[div^.w] dx+ / $^(1 - G(-A£,7v))[div^.w] dx, 
where, as stated in (16.11) . 

|t ^ <l>eG'(-A^,jv)[div^.w] dxj ^ in L'^{0,T). (6.6) 
Writing $^ = $i + $2 as in ( K23\i . (gJl]), we have 

while, in agreement with (I5.17P and Lemma 15. ![ 

|t ^ $^(1 - G'(-A,,^))[div,w] da;| (6.8) 

-> |th^^$^(l-G(-AAr))[div^w] dxj inL2(0,r), 
where the resulting expression is small as soon as G ~ l[o,oo)- Indeed 

f $1(1 _ G(-A^))[div,w] dx= [ (-A^.)i/2$i— 1--(1 - G(-A;v))[div,w] dx, 
Jn Jn {—Am) ' 



where 
while 







(-A.)V2t^"^"^] 


L2(n) 



AA.)^/^[$^]|U.(n) = ||V..<f^|U.(f.), 

< ||w||L2(n;_R3). 



L2{n) 



V,(-Aiv)"'[div,w] 



L2(f7) 



Relations (16.31 - [6^ imply (16. 2p . in particular, we have shown (I2.23p . The proof of 
Theorem 12.11 is now complete. 
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6.2 Boundary behavior of the Umit velocity field U 



In the previous analysis, we left open the problem of the boundary conditions satisfied by 
the limit velocity field U. To this end, we revoke the results of |8]. Suppose that, after 
suitable translation and rotation of the coordinate system, a part F of the boundary of Q 
can be described by a graph of a function b G W^'°°{U), U C R^, 

r = {(xi, X2, X3) I (xi, X2) eU, xs = b{xi, X2)}, 

while = dfls riU x R are represented as 

= {{xi,X2,Xs) I (xi,X2) eU, X3 = b^{xi,X2)}, 

where {&e}£>o is a bounded sequence in W^'°°{U), 6^ — > 6 in C{U). Similarly to [8], we 
assume that the boundaries are oscillating for e — 0. More specifically, introducing a 
Young measure 7l[y], y & U, associated to the gradients {'Vybe}e>o, we suppose that 

supp[7^[?/]] contains two independent vectors in R^ for a.a. y E U. (6.9) 

As shown in [S], condition ( 16 ■9p implies 

U|r = 0, 

see also Bfezina [5] for refined results in this direction. 
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